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The large iV limit of the hermitian matrix model in three and four Euclidean space-time 
dimensions is studied with the help of the approximate Renormalization Group recursion formula. 
The planar graphs contributing to wave function, mass and coupling constant renormalization 
are identified and summed in this approximation. In four dimensions the model fails to have 
an interacting continuum limit, but in three dimensions there is a non trivial fixed point for the 
approximate RG relations. The critical exponents of the three dimensional model at this fixed 
point are v = 0.665069 and r\ = 0.19882. The existence (or non existence) of the fixed point and 
the critical exponents display a fairly high degree of universality since they do not seem to depend 
on the specific (non universal) assumptions made in the approximation. 
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1) Introduction. 



Since 't Hooft's original proposal [1], the large N limit has been regarded as one of the most 
promising attempts to understand strongly coupled gauge theories. The observation that a field 
theory simplifies in the presence of a large symmetry group is, in fact, more general and has led, 
among other things, to an understanding of many linear and non linear a models in dimension 1 ) 
2 < d < 4 [2] , to the exact solution of matrix models in dimension d < 1 [3] and to the solution of 
d = 2 QCD [4]. 

The great simplification arising in the study of a models is the vector-like nature of the fields, 
i.e. the fact that the number of degrees of freedom grows like N. This allows one to reduce the sum 
of all the leading "cactus" diagrams to the sum of a finite number of effective graphs, providing in 
this way a tractable perturbative expansion. The leading term in the expansion already captures 
the relevant features of the various models, including spontaneous symmetry breaking (or the lack 
of it in d = 2), confinement, dynamical generation of relevant operators, etc... 

The situation is more complicated in the case of matrix models, of which QCD is the prime 
example. There, the number of degrees of freedom grows like N 2 and one is faced with the arduous 
task of summing the leading "planar" diagrams. This problem has been solved in [5] for the single 
matrix model (d = 0), and for the quantum mechanical anharmonic oscillator (d = 1), with its 
remarkable equivalence to a free Fermi gas. These low dimensional models have been shown to 
describe d = 2 quantum gravity and non critical string theory [3]. 

Alas, progress has been much slower in the study of higher dimensional matrix models. The 
factorization properties of the expectation values of observables [6] has led to the conjecture [7] that 
the leading order in 1/N should be dominated by a single field configuration, the "master field". 
Very recently the precise meaning of the master field has been clarified using non commutative 
probability theory [8]. This recent development is very exciting but is also somewhat worrisome: 
the connection between the master field and the knowledge of the connected Green functions is 
so explicit that it might signify that finding the master field is just as hard as exactly solving the 
large N theory without it! It seems likely that one would have to develop approximations to this 
construction, treating the large N theory as an "exact" theory in its own right. 

In this paper, we develop an approximation to d = 3 and d = 4 hermitian matrix models with 
quartic potential using one of the oldest tools in the non perturbative formulation of quantum 
field theory: the approximate renormalization group (RG) recursion formula [9]. We make some 

Throughout this paper, the letter d refers to the dimension of (Euclidean) space-time. 
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modification to this technique to make it more suitable for the problem at hand: on the one 
hand, we relax the restriction of no wave function renormalization [10] to allow for a non trivial 
anomalous dimension ij and, on the other hand, we restrict our attention to the renormalization 
of the quartic term in the interaction. The first change is crucial because, to leading order in l/N, 
there are already an infinite number of diagrams contributing to 77; neglecting them would ruin the 
structure of the RG flow. The second assumption is not crucial, but it simplifies the analysis by 
allowing one to perform explicitly the integrals in the recursion formula using the results from the 
single d = matrix model and to reduce the single recursion step from an integral equation to an 
algebraic one. 

The main results of this investigation are as follows. In d = 3 the model has a non trivial fixed 
point in addition to the Gaussian one. The large TV critical exponents at this point are computed 
to be r\ = 0.19882 and v = 0.66517. In d = 4 the model fails to have a non Gaussian large ./V limit; 
its continuum limit is a free theory. This last result may very well be due to the roughness of the 
approximation and should not discourage us to study the d = 4 theory further. 

The paper is organized as follows. In section two, we present the model and qualitatively de- 
scribe its renormalization in the large N limit. In section three, we give a sketch of a perturbative 
calculation to third order in the coupling constant. This section is not needed for the further calcu- 
lations which are intrinsically non-perturbative and it is only given for comparison. In section four, 
we derive the approximate recursion formula. All the assumptions that go into its derivation are 
spelled out and explicit expressions for wave function, mass and coupling constant renormalization 
are obtained in the large N limit. In section five, we study the RG equation derived, establish the 
existence of a non Gaussian fixed point for d = 3 and compute the critical exponents. We also 
comment on the non existence of a non Gaussian fixed point in d = 4. All the results for the d = 
hermitian model [5] of direct relevance for this paper are collected in the appendix. 



The model we study in this paper is the d = 3, 4 hermitian matrix model with quartic 
interaction. The field variable is a 7V x N hermitian matrix $(#) = In the spirit of Wilson's 

RG, we assume the presence of an effective cut-off A in momentum space. Choosing the mass 
unit to be A itself, allows us to set A = 1 and represent all quantities in dimensionless units. In 
this way we write the Fourier transform of the field as 



2) The hermitian matrix model and its RG. 




(1) 
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The Euclidean action for reads 



S[$\ = ^Jd d x tr (^(W) 2 + r $ 2 ^j + ^ f d d x tr <S> 4 

u /" d d k 1 d d k 2 d d k 3 / ., > . , , , N 

^ io<ifc u , 3 |<i ( 2 ^r V 



'0<|fcl,2,3l 

It is well known that, in the l/N expansion, the leading term for the free energy is of order 
0{N 2 ). Also, the correct normalization for singlet operators, in order for them to have a finite 
vacuum expectation value, is as follows: a factor of 1/y/N for every power of the field and a factor 
of l/N for every trace; e.g.: 

1 tr$ fc - tr <£ fc tr ••• (3) 



jyl + fc/2 ' jy 2 + (fc+Z)/2 

Action (2) is, of course, just a truncation of the most general action in theory space consistent 
with the symmetries of the problem (throughout the paper we deal with a parity preserving theory) . 
A generic term in the most general action would contain products of traces of powers of the field 
multiplied by (short ranged) functions of their momenta and by the appropriate powers of l/N: 



1 1^ " ' U{h "' )tT • • • ) tr •••)•••• ( 4 ) 



1 f d d h 



That (2) is a consistent truncation can only be checked at the end. However, it should at 
least seem reasonable because it includes the three most relevant operators at the Gaussian point. 
Actually, there is a fourth such operator: (v /N 2 ) [ d d x[ tr $ 2 ) 2 . One might be tempted simply 
to dismiss such an operator as trivial because it can be removed by introducing an auxiliary 
singlet field A to the action —(N 2 /v ) J d d x[\ — (v /N 2 ) tr ($ 2 )) 2 , or because it only affects 
the correlation functions of the theory through the insertion of vacuum diagrams via a contact 
interaction. However, this is not quite enough; if such a term was generated by the RG it would 
affect mass renormalization and its effect would have to be taken into account. 

At the cost of being pedantic, we begin by showing that such a term in not generated by 
applying the RG transformation on (2) and therefore it may be excluded a priori. This boils down 
to showing that actions of the form (2) are mapped into themselves by the RG, modulo terms of 
higher Gaussian dimension. There might be some confusion on this point since it is well known that 
the usual exponential relations between the generating function of the connected and disconnected 
Green functions, as well as the Legendre transformation between the generating function of the 



connected Green function and the one particle irreducible vertices, are not valid in the large N 
limit. Why then should the Wilsonian action re-exponentiate on itself without generating terms 
non linear in the trace? The solution to this puzzle is that there are leading Green functions in 
1/N that are not planar (fig. 1) and, while they do not appear in the generating function, they do 
contribute to the Wilsonian action. 




+ 



+ 



X Q O 




Fig 1.) The usual decomposition of the four point function into its connected components. The 
last contribution is not "planar" but should be included when computing the renormalized action. 



Let us address this and other related issues by carefully carrying on the first few steps in 
the renormalization programme. We must integrate out the higher momenta from the action and 
obtain a new action that, when rescaled, will yield information about the RG flow. Let us then 
begin by splitting the field variable into slow and fast components $ = $ + with momenta smaller 
and larger than 1/2 respectively: 



$(x) = f 



d d k 

'0<|fc|<l/2 ( 27r ) d 

Action (2) decomposes into 



e ix - k $(k) and 



m = [ 



d d k 

l/2<|fc|<l ( 2 ^) d 



ix-k 



4>{k)- 



(5) 



S[$]=Sm+a[$,cf)}+S[cj>}, (6) 
where S[$>] and S[4>] are the same as (2) with the appropriate restriction on the momenta and 

U ° ' l'l i " / I ,,T V S , I. 2,^2 , oi. ,U;2 , ,,,i.r. i (J) 



a[&, = ^ J d d x tr + 4^$ 2 + 2{(j)$>y + # 3 $J . 

(Of course, the quadratic term factorizes by momentum conservation). Let us denote, for any local 
functional <j>] of the fields, 

fV^e-^F 
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(8) 



The denominator in (8) is such that < 1 >= 1. It removes all vacuum diagrams just as in the 
ordinary case. The relation between Green functions with and without vacuum diagrams is the 
familiar one because every leading vacuum diagram is planar and vice versa. 

Expanding exp(— a) to forth order in <£> and keeping only those terms that are allowed by 
parity one obtains 

^e"^ =^1-^ y^2^4tr0 2 $ 2 (x) + 2tr (0) 2 (x)^j 

+ \w j ddxdd y{^ 2 tr <l& 3 ( x ) tr <t?Hy) + 16 tr 3 $(x) tr 3 $(y) 

+ 16 tr <f) 2 $ 2 (x) tr <f) 2 ® 2 (y) + 16 tr 2 $ 2 (x) tr (0) 2 {y) + 4 tr (#) 2 (x) tr (0$) 2 (y)^ 

- j d d xd d yd d z ^192 tr 2 <l 2 (x) tr <j?${y) tr 3 $(z) 

+ 96 tr (#) 2 (x) tr 3 <%) tr 3 $(z)^) 

+ 2lmJ ddxdd V ddzddw {^^ tr <t?Hx) tr 3 ^(2/) tr 3 <l(z) tr 3 ^(w)^j ^. 

(9) 

(This is not a perturbative expansion; each term contains diagrams with an arbitrary number of 
(f) 4 vertices.) 

The first thing to notice is that none of the four terms containing the operator tr {4>§) 2 in 
(9) contributes to the expression. Terms like those in fig. 2. a yield a contribution proportional 
to (1/N)( tr l>) 2 . This term would be leading (0(N 2 )) if tr l> was allowed to pick up a vacuum 
expectation value, but it is subleading in a parity invariant theory: 
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1( tr *) 2 ) = N 2 x ^( tr $) 2 ) = N 2 x ( (JL. tr ^ + O^" 2 ) 
N 2 x (0 + O(iV- 2 ) ) =0(1) 



(10) 



Terms like those in fig. 2.b "would be" leading but they vanish identically because of momentum 
conservation. In fact, all diagrams containing a sub-diagram with one i> leg and one 4> leg (fig. 3) 
vanish identically by momentum conservation. This implies that the term tr <ft<b 3 (x) tr ^> 3 <i(y) in 
(9) (fig. 2.c) also vanishes identically. Finally, notice that a term like the one in fig. 2.d is planar 
but it is subleading because the fields $ are not connected at the outside of the diagram. 




Fig 2.) Various contributions to the renormalized action in t' Hooft's double line notation. All 
external lines represent the slow field $ and all the internal lines the fast field <ft. 2. a) is subleading 
if the trace does not pick up a vacuum expectation value. 2.b) and 2.c) are identically zero by 
momentum conservation. 2.d) is planar but subleading. 




Fig 3.) Same notation as in fig 2. All diagrams containing this sub-diagram vanish. 
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^ I d d x( tr</> 2 <l 2 

32tt / ,7'/ „ 7'/. ..;</. / , .. ..2^2, ,.\ ,,. .Ji,Tw ,» , .. 

iV 3 



3^0/ / „/...,/ 
32u 4 



Considering all the remaining terms and performing the decomposition into connected com- 
ponents one obtains: 

/e-M =1-^1 d d xl tr </> 2 <l 2 (x)\ + ^§ J d d xd d yl tr </> 3 $(x) tr ^(y)\ 

\ I \ I conn. \ / conn. 

+ ^ C / d d xl tr 2 l> 2 (x)\ ) + ^ / W tr ^ 2 $ 2 (x) tr 2 l> 2 ( y )\ 

V \ / conn./ \ / conn. 

(x) \ J d d yd d zl tr c/> 3 <l(y) tr ^ 3 $(z) \ 

/ conn. \ / conn. 

d d xd d yd d z( tr 2 $ 2 (x) tr 3 <%) tr 3 $(z) \ 

\ / conn. 

J d d xd d yl tr 3 $(x) tr 3 $(y) \ J 

\ / conn./ 

y d d xd d yd d zd d wl tr 3 $(x) tr (/> 3 $(y) tr (j) 3 $(z) tr </> 3 $(t(;)\ 

\ / conn. 

(11) 

The only relevant point is that the next to last quantity ( f d d xd d y( tr <p 3 ^(x)(f> 3 ^(y)) conn ) 2 is 
multiplied by the factor 32 = 3 x 256/24, in other words we had to use all three disconnected Green 
functions in fig. 1 since we were taking the vacuum expectation value of four singlets ( tr </> 3 l>). 
This makes it possible to re-exponentiate (11) without generating terms non linear in the trace. If 
we denote by <!?[$] the exponential: 

e -sm = ^ e -^y (12) 

then 

5[$] =^ J d d xl tr 2 $ 2 (x) \ - M y tr 3 $(x) tr </> 3 <%) \ 

\ / conn. \ / conn. 

y d^y/tr^a) tr^* 2 (y)\ 

\ /conn. (13) 

^ y d d xd d yd d zl tr 2 $ 2 (x) tr 3 <l(y) tr 3 $(z) \ 
32u 4 



3iV 4 



£ y d d xd d yd d zd d wl tr 3 $(x) tr 3 $(y) tr </> 3 $(z) tr 3 $H 



So far, our discussion has been completely general and no approximation has been made, 
except for the truncation to operators containing no more than four fields Of course, the 
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identification and evaluation of the diagrams contributing to (13) is still a formidable problem and 
we must make further approximations; these will be made in section four. Before that, in the next 
section, we present the sketch of a perturbative calculation as a warm up. 

3) A perturbative calculation. 

A perturbative expansion to (0(uq)) around the Gaussian point and to leading order in 1/N 
results in the diagrammatic expansion for given in fig. 4. 




Fig 4.) Perturbative expansion of the renormalized action in single line notation. The external 
lines represent the field <& and the internal ones the field <p. 
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All diagrams that are forbidden by momentum conservation have been omitted. The symmetry 
factors are all powers of 2 and this is typical of planar expansions. One can compare these factors 
with those in [5], keeping in mind that there is an overall factor 2 of difference in the two-point 
function and an overall factor 4 in the four-point function. This is so because, here, the external 
legs are attached to the slow field l> acting as a background field. One must also remember the 
explicit factor of 1/N in front of the four-point functions. 

Just as an example, the first diagram in fig. 4. corresponds to the expression 

CM)"^ <"> 

whereas, say, the tenth diagram corresponds to 

_M f d d xd d y tr (&(x)$ 2 (y))G(x,y) 2 . (15) 

In (14) and (15), the quantity G(x,y) is given by the tree level propagator for (f>: 

G(x,y) =(n~ 2 tr cj>(x)cj>(y)) =[ ^L^-v) G {k) 
\ io Jo<\k\<i (2v!T 

' for \k\ < \ (16) 

G(k) = < 

[ WTF for | < < 1 

The symbol ^ ■ ■ ■ ^ is the Gaussian average and one can also write, for each component of the 
matrix, 

^(aO#(y)^ = S?qG(x, V ) + 0(1). (17) 

One could then perform the integrals and obtain a perturbative expression for S, but this 
will not be done explicitly because the main purpose of this section was to present the relevant 
diagrams with their relative weights for later comparison. Let us therefore go on to the next section 
and present those approximations that allow one to obtain non perturbative results. 

4) The approximate RG formula in the large N limit. 

We must now make some different approximations in order to compute (13) without using 
perturbation theory. The problem is, of course, that in the expression for S we want to keep the 
contribution of integrals of arbitrary higher number of loops and we do not know how to do it 
exactly. The first thing we do is to set all incoming momenta (the momenta of the slow field <l) to 
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zero, in order to obtain the "ultra local" contributions to S. This approximation is equivalent to 



the truncation we have chosen for (2) and it is not an independent assumption, but it should only 
be applied to those terms that are already ultra local in (2), namely, mass term and interaction. If 
we applied it to the kinetic term we would rule out the possibility of anomalous dimension for the 
field. This would be too strong since, in the large N limit, there is an infinite number of planar 
diagrams contributing to wave function renormalization. Let us therefore start with the mass and 
the interaction terms, leaving the analysis of the kinetic term for later. 

In the spirit of the approximate recursion formula, we now make two further assumptions: 

1. Replace every propagator (regardless of its momentum dependence) by 1/(1 + r ) 

2. JReplace every loop integral by a constant Cd 

The actual value of the constant Cd is not important. To be specific, we will assume 



Now we must identify those diagrams that contribute to S and estimate the relative powers 
of Uq, 1/(1 + ro) and Cd in each of them. Since all external momenta are set to zero, momentum 
conservation implies that not only diagrams containing sub-diagrams like the one in fig. 3, but 
all diagrams that can be disconnected by cutting one internal 4> line vanish, no matter how many 
external $ lines are attached to it. Therefore, only one particle irreducible diagrams contribute, as 
one would naively have expected. 2 ) The great simplification in considering only planar diagrams 
is that, with the assumptions above, the sum of these diagrams can be explicitly evaluated by using 
the single matrix model (c.f.r. appendix). None of the further assumptions usually made in the 
approximate recursion formula, (such as forbidding those diagrams with an odd number of internal 
lines at some vertex,) need to be made. The situation is summarized in fig. 5. The diagram in fig. 
5. a is exactly forbidden; the diagram in fig 5.b is forbidden under the assumption of zero incoming 
momentum; the diagram in fig. 5.c, that would be forbidden by the ordinary approximation, is in 
fact allowed, as it should. 

2 ) This observation is not directly relevant to the calculation of the corrections to r and u 
because, due to the particular nature of quartic interactions, all one particle reducible diagrams 
with two or four external legs have a sub-diagram of the kind in fig. 3; but it should be kept in 
mind if one wants to compute diagrams with more than four external legs. 




0.014776 for d = 3 
0.002968 for d = 4 ' 



(18) 
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Fig 5.) Some planar diagrams in single line notation. 5. a) has already been shown to be identically 
zero. 5.b) vanishes in the approximations used in this section. 5.c) does not vanish and it is 
accounted for in our formulas. 



To estimate the relative powers of u , l/(l+r ) and ca, we make use of the following topological 
formulas: for a particular diagram with V vertices, E > external l> legs, P propagators and L 
momentum loops, the following two relations apply 

AV = E + 2P and V — P + L = 1. (19) 

The reader should not be confused by the second relation in (19). In (19), L refers to the number 
of loop momentum integrals, not color loops. Among all diagrams with L momentum loops, the 
planar limit picks out those with the maximum number L™^ x of color loops. For E > these two 
numbers coincide, whereas for vacuum diagrams i™^ x = L + 1. 

For a given number E of external legs, we can solve for P and L as function of the number of 
vertices V. We can therefore write the mass term (E = 2) S ro and the interaction term {E = 4) 
S UQ in S as follows 

S ro W = (l + ro)fJj^^] fd d x tr l> 2 (x) 

^} = pu (jT^J jd d x tr ^( X ), 

where f ro (g) and f Uo (g) are given in terms of the single matrix model proper vertices (c.f.r. ap- 
pendix) : 

sm = T -^r 1 «i /«,(,) = T -Mz±. (2i) 
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After all that has been said, the proof follows easily. First, we match the powers of g in the single 
matrix model with the powers of u in the full quantum theory by dividing by the overall factors 
of 2 and 4 already encountered in perturbation theory. 3 ) Second, we subtract the tree level terms, 
already present in (6) and not to be included in S[&] (13). Finally, we match the powers of 
Cd and 1/(1 + r ) with the help of the topological formulas (19). 

We must now face the problem of wave function renormalization by relaxing the condition of 
zero external momenta. Let q « 1/2 be the external momentum of the fields First we isolate 
those diagrams that contribute to the renormalization: they are those with two external legs 
connected at two different points in space-time. Terms where the two external legs are connected 
to the same point do not have any explicit dependence on the external momenta and they should 
not be included. The decomposition of the one particle irreducible two point function into these 
two kind of terms is given in fig. 6. 




Fig 6.) A useful decomposition of the vertex T2- Solid blobs represent connected Green functions 
and crossed blobs represent one particle irreducible vertices. The last term contributes to wave 
function renormalization. 

Fig. 6 can be checked explicitly for the single matrix model (c.f.r. appendix): 

T 2 (g) = 1 + 8gG 2 (g) - AgT 4 (g)G 3 2 (g). (22) 

Only the last term in (22) contributes to wave function renormalization. To find such contribution 
we must take the derivative with respect to q 2 of each contributing diagram and then set q 2 = 0. 

3 ) The explicit factor of g at the denominator in the definition of f Uo is there simply because we 
have factored out one power of u in the definition of S uo , it does not represent a singularity at 
the origin! 
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Again, this is hard to do exactly because there may be more than one q dependent propagator 
in a diagram. We overcome this difficulty by making the further approximation, (justified by 
dimensional analysis and also by acting with the derivative on a single propagator): 



_d 

dq 2 1 + r 



■ . ( 23 ) 



This allows us to estimate the planar contributions without any effort. It is an approximation in 
the same spirit as those made before. It should correctly capture the qualitative behavior and the 
order of magnitude of the correction. By using it, we obtain the last piece of S: 

§ A*\ = U ((TT^ji) / ddx tr ( v< ^)) 2 > ( 24 ) 

where (remembering the 1/2 normalization factor as before), 

f q2 (g)=2gT 4 (g)G 3 2 (g). (25) 
Action (13) is therefore specified in terms of three known functions f ro , f Uo and j q i of 

These three functions can be written as rational functions of a single function of g (c.f.r. appendix) 

A S ) = ^^zl (27, 



as follows 



, N (a 2 - l)(a 2 -3) o o 

fr ° i9) 2a 2 (4-a 2 ) * i9 ~ 4 ° 9 + 8329 + 
27(5 -2a 2 ^ 
~ (4 -a 2 ) 4 
(q 2 -l) 2 (5-2a 2 ) ^ ^ 2 or/ . , 
18a 2 (4-a 2 ^ 



^ (9) = ^'^_^y - 1 ~ ~ 8 9 + 22 V - 7296/ + • • • (28) 

, , (a 2 -l) 2 (5-2a 2 ) , , 
fA9) = ] « 8g 2 - 256/ + • • • . 



A plot of these three functions is given in fig. 7. The range in which they are real is— 1/48<<7<oo; 
they are analytic at g = and their expansion given in (28) matches the perturbative expansion 
in fig. 4 if one makes the the same approximations as we have made here. 
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Fig 7.) A plot of the three functions described in the text. f ro is the mass renormalization, f uo 
the coupling constant renormalization and f q 2 the wave function renormalization. f Uo approaches 
a finite limit as g — > oo whereas f ro and f q 2 diverge as y/g. 



Let us summarize what has been done so far. We started with the action S[$>] (2), and after 
integrating out the fast modes 4> we obtained a new action for the slow fields <£: 

2 



= + = ^jd d x 



1 + 2/, 



tr V$(x) 



+ - / d a x 



r + 2(l + r o)f J^ 



tr + jj d d x 



1 + / 



«0 



(l + r ) 2 



tr $(x) 4 . 
(29) 

To obtain the RG equations we have to rescale co-ordinates x — ► 2x, (A; — > fc/2) and redefine 
a new field 

*(*) = 2 (d " 2)/ yi + 2/ g2 ( IT ^)$(2x) (30) 

in terms of which the action becomes again of the form (2), but with a new mass term and a new 
coupling constant 

S[§\ = ^Jd d x tr (V$) 2 + n$ 2 ^ + ^ J d d x tr $ 4 . (31) 
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This is just the first step of the RG transformation. The same relation between r ,u and ri,U\ 
holds for a generic step between r;, ui and r; + i, ui +1 . We have therefore obtained the RG equations 
for the hermitian matrix model in this approximation: 

n + 2(1 + n )f ro (g) 
i + VAa) 

where g = g(ui,ri) is given in (26). In the next and final section we investigate these equations to 
determine the dynamics of the theory. 

5) Study of the RG. Fixed point and critical exponents. 

The iteration of eqs. (32) determines the RG flow in the two dimensional space spanned by 
two co-ordinates r and u. At a fixed point r*,u*, any further iteration leaves the co-ordinates 
invariant, i.e., r* = r% = r/ + i and u* = ui = ui + \. Let us also define, for convenience, the quantity 

9* = 9(u*,r*) = j^y 2 . (33) 



At the fixed point, eqs. (32) then become 

r* =4 



r* + 2(1 + r *)f ro (g*) 
l + 2/ 2 (g*) 



,* _o4-cL.* l + /«oU?*) 



u* =2 4 "<V 



(l + 2f q2 (g*)y 

By eliminating the solution u* = (Gaussian point) from the second eq. in (34), we obtain an 
algebraic equation is terms of the variable g* alone: 

2 d -\l + 2f e (g*)Y = l + f UQ {g*). (35) 

For d = 4, one can see right away that the only solution is g* = 0. In fact, for g > (physical 
region), it is always f q 2(g) > and f Uo {g) < (fig. 7), thus preventing (35) from having non trivial 
solutions 4 ) . Within our approximation, this lack of a non trivial fixed point in four dimensions is 
rather robust since it only depends on the relative sign of f q 2 and f Uo . However, it should not be 
interpreted as a "no go" result because more subtle approximations might reveal a richer structure. 

Things are much better in d = 3. We can immediately see that eq. (35) must have a non 
trivial solution in the physical region by considering the asymptotic behaviour of the two sides of 



4 ) One can also check that there are no solutions in the interval [— 1/48, 0[ 
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the equation as g — > and g — > oo. As g — > 0, the l.h.s. — > 1/2 and the r.h.s. — > 1, whereas 
as g — ► oo, the r.h.s. approaches a finite value and the l.h.s. blows up. The two functions must 
then cross somewhere. Numerically, the only two solution in the whole interval [—1/48, +oo[ are 
the Gaussian point and g* « 0.408655. By plugging in this last value to the first of (34) one 
obtains r* = —0.642902 and, recalling the definitions of g* (33) and of ca (18) one also obtains 
u* = 3.52675. 

If one studies the evolution of this fixed point as d goes from 3 to 4, one finds that it approaches 
the Gaussian one. In this sense, our fixed point is " Wilson-Fisher- like" . On the other hand, our 
result does not require using the e expansion, and it is a truly large N result because it hinges on 
the sum of planar graphs and the consequent asymptotic behavior of the functions f ro , f Uo and f q 2 
for large g. As a final comment, notice that, if we neglect wave function renormalization by setting 
f q 2 = 0, the fixed point goes away, yet another indication of the different nature of the calculation. 

To compute the large N critical exponent v we linearize the RG equations (32) near the fixed 
point. Numerically: 

r i+1 - r* =1.97107(n - r*) + 0.303152(w z - u*) 

(36) 

u l+ i - u* =6.11699(n - r*) + 0.690315(u z - «*). 
The largest eigenvalue A max . of the matrix 



/ 1.97107 0.303152 \ ( 
^6.11699 0.690315/ [ ' 

= 0.665069. (38) 



M 

is A max . = 2.8355, yielding a critical exponent 

log 2 
v = 

log A m£ _ 

As a " check" of universality, one can compute the dependence of M on the non universal quantity 
Cd and notice that, while the off diagonal entries of M depend on Cd, the characteristic polynomial 
(i.e., trace and determinant) does not; and therefore neither does v. 

To compute the anomalous dimension r\ of the matrix field we simply express the wave function 
renormalization as a power of 2 and take the logarithm: 

log(l + 2/ 2(o*)) , « 

= s v t j q \y i) = o_i9882. 39 

log2 V ' 

The exponent r\ is also independent on cj. 
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Appendix: Some useful facts about the single hermitian matrix model. 

In this appendix we collect, without proof, some basic facts about the single matrix model 
that are used throughout the paper. All details can be found in [5]. 

Let 4> be an N x N hermitian matrix. The action of the d = hermitian matrix model we are 
interested in is 

S d =°[<f>]= tr (Itf + lA. (40) 



For every function F{4>) we define 

v d=0 

' AT— kx> fd N2 <l>e- sd=0 

The Green functions are defined as 

d=0 



m) =hmL-^— — (41) 



^(2) = (^tr^ + 0(±), (42) 

where the factors of 1/N are chosen so that the leading (planar) contribution is finite. 

The relation between the two and four point Green functions G 2 and G4, the connected ones 
C 2 and C4 and the vertices T 2 and T 4 is 

C 2 = G 2 

C4 = G4 — 2G 2 



(43) 

2 — 2 — Lr 2 



To = Co 1 = G 1 



Y4 — — G±C 2 ^ — (2G2 — G^)G 2 ^ 
The 5 dependence of all planar Green functions and vertex functions can be written in terms 
of rational functions of a 2 (g) where 



.4,2 1 n :„ „2, „x VI + 48ff - 1 



12 5 a 4 + a^ - 1 = i.e., a\g) = v ^ . (44) 



In particular, for the Green functions of interest: 



G 2 = -a 2 (4 - a 2 ) « 1 - 85 + 144c/ 2 - 3456# 3 + • • • 
G 4 = a 4 (3 - a 2 ) « 2 - 36 5 + 864c/ 2 - 24192a 3 + ■ ■ ■ 
C 2 = \a 2 {A - a 2 ) « 1 - 85 + 144c/ 2 - 3456# 3 + • • • 
a 4 

C 4 = - y (1 - a 2 )(5 - 2a 2 ) « -4a + 160# 2 - 5760 5 3 + 

r 2 = 2TT- 2 ^ ~ 1 + 8<7 - 80 5 2 + 1664<7 3 + ■ ■ ■ 
a^(4 — a 2 ) 



(45) 



9( l-a 2 )(5- 2a 2 ) 
a 4 (4-a 2 



r 4 = 2U ~ 4 5 - 32 5 2 + 896 5 3 + 
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